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Abstract 

Transmission through and reflection from a potential barrier, and the very closely 
related issue of particle production from a parametric resonance, are topics of consid- 
erable general interest in quantum physics. We have developed a rather general bound 
on quantum transmission probabilities, and recently applied it to bounding the grey- 
body factors of a Schwarzschild black hole. In this current paper, we take a different 
tack - we report a way of using the Miller-Good transformation (which maps an initial 
equation to a final equation for a different potential) to significantly generalize the 
previous bound. We then apply this general formalism in a very specific manner to 
derive a rigorous bound that is "as close as possible" to the usual WKB estimate for 
barrier penetration. 
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1 Introduction 

Following the discussion presented in [TJ, consider the Schrodinger equation, 

u(x)" + k(xfu(x) = 0, (1) 

where, k(x) 2 = 2m[E — V(x)]/h 2 . As long as V(x) tends to finite (possibly different) con- 
stants V±oo at left and right spatial infinity, then for E > max{V +00 , V-oo} one can set up 
a one-dimensional scattering problem in a completely straightforward manner - see for ex- 
ample any of the standard references [2l[3llll[5l[6l>[7|,[8],[9j US] - The scattering problem is 
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completely characterized by the transmission and reflection amplitudes (t and r), though the 
most important aspects of the physics can be extracted from the transmission and reflection 
probabilities (T = |t| 2 and R = \r\ 2 ). Relatively little work has gone into providing general 
analytic bounds on the transmission probabilities, (as opposed to approximate estimates), 
and prior to the work leading to the current article, the only known result as far as we have 
been able to determine is this [UJ [T2] : 

Theorem 1. Consider the Schrddinger equation (QJ). Let h(x) > be some positive but oth- 
erwise arbitrary once - differentiable function. Then the transmission probability is bounded 
from below by 

( ,co J (h') 2 + (k 2 - h 2 ) 2 ) 

^"{jL^^-is ~ d j- (2 » 

To obtain useful information one should choose asymptotic conditions on the function h(x) so 
that the integral converges - otherwise one obtains the true but trivial result T > sech 2 oo = 0. 
(There is of course a related bound in the reflection probability, R, and if one works with 
the formally equivalent problem of parametric oscillations, a related bound on the resulting 
Bogoliubov coefficients and particle production.) 

This quite remarkable bound was first derived in [TT] . with further discussion and an 
alternate proof being provided in [12]. These bounds were originally used as a technical step 
when studying a specific model for sonoluminescence [13] , and since then have also been used 
to place limits on particle production in analogue spacetimes [H] and resonant cavities [T5] . 
to investigate qubit master equations [16], and to motivate further general investigations of 
one- dimensional scattering theory [17]. Most recently, these bounds have also been applied 
to the greybody factors of a Schwarzschild black hole |18j . 

A slightly weaker, but much more tractable, form of the bound can be obtained by 
applying the triangle inequality. For h(x) > 0: 



T > sech 2 



ln(/i)'| + 



\k 2 - h 2i 
h 



dx \ . (3) 



Five important special cases are: 

1. If we take h = koo, where koo = lim^^ioo k(x), then we have [TJ [Til H2 



r-OO 

T>sech 2 <i^/ \kl~k 2 \dx). (1) 



"oo J — oo 



2. If we define k +OQ = \im x _> +00 k (x) ^ k^^ = lim x ^_ 00 k(x), and take h(x) to be any 
function that smoothly and monotonically interpolates between /c_oo and k +OQ , then we have 



T > sech 2 



In 



&+oo 



1 r°o I*; 2 -/* 2 ! 
+ o / J 1 - dx \ ■ ( 5 ) 



— oo 
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This is already more general than the most closely related result presented in [TT1 [T2] . 
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3. If we have a single extremum in h(x) then 



T > sech^ 



In 



oo' v — oo 
lb ext 



1 Z" 00 \k 2 -h 2 \ , . 

^ . (6) 



2 ./_„ /i 



This is already more general than the most closely related result presented in [TTI [T2] . 
4. If we have a single minimum in /c 2 (x), and choose h 2 = max{fc 2 , A 2 }, assuming fc^ in < 
A 2 < fc^, (but still permitting fc 2 iri < 0, so we are allowing for the possibility of a classically 
forbidden region), then 

T > sech 2 (1 In ( t±22^) + _L [ | A 2 - fc^l . (7) 



2A 



A 2 >fc 2 



This is already more general than the most closely related result presented in [TTI [T2] . 

We emphasize the A is an adjustable parameter that can in principle be used to optimize 
the bound - furthermore it is this special case that has a close connection to the WKB-like 
bound we shall eventually derive later on in this article. 
5. If k 2 (x) has a single minimum and < fc ^ iri < k 2 .^, then 

T>sech 2 |-ln^ +ooA; - 00 



k 



2 

min 



This is the limit of (J7|) above as A — > k m \ n > 0, and is one of the special cases considered in 

Comment: In the current article we shall not be seeking to directly apply the general bound 
02]), its weakened form <^j, or any of its specializations as given in (SJ)- (JHJ) above. Instead we 
shall be seeking to extend and generalize the bound to make it more powerful. 
The tool we shall use to do this is the Miller-Good transformation [TJ IT9] . 



2 The Miller-Good transformation 

Consider the Schrodinger equation (jTJ, and consider the substitution [TJ 

u{x) = -yl=U(X(x)). (9) 

We will want X to be our "new" position variable, so X(x) has to be an invertible function, 
which implies (via, for instance, the inverse function theorem) that we need dX/dx ^ 0. In 
fact, since it is convenient to arrange things so that the variables X and x both agree as to 
which direction is left or right, we can without loss of generality assert dX/dx > 0, whence 
also dx/dX > 0. 

Now compute (using the notation Ux = dU/dX): 

u\x) = U X (X)VX- l --^—U{X), (10) 
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and 

I Y"> q (X") 2 

u"{x) = U xx {X){X'f 12 - - 2jYW2 U + ~ 4 ^j- 2 U. (11) 

Insert this into the original Schrodinger equation, u(x)" + k(x) 2 u(x) = 0, to see that 

( k 2 1 X'" 3 (X") 2 1 

fa+ W"w + Jwr =0, (12) 



which we can write as 
with 



U xx + ^C/ = 0, (13) 

K* = J-i k *-l^ + 3 ~^\ (14) 

(X') 2 1 2I' + 4 (X') 2 J ' 1 J 

That is, a Schrodinger equation in terms of u(x) and fc(x) has been transformed into a 
completely equivalent Schrodinger equation in terms of U(X) and K(X). You can also 
rewrite this as 

K* = 7 ^lh? + >&(^=Y\. (15) 



(X') 2 I \Vx> 

The particular combination 

shows up in numerous a priori unrelated branches of physics and is sometimes referred to as 
the "Schwarzian derivative" [20jJ3 

• As previously mentioned, to make sure the coordinate transformation x -<->• X is well 
defined we must have X\x) > 0. Let us call this j(x) = X'(x) with j(x) > 0. We can 
then write 

l/ + 3(/) 2 
f I 2 j 4 f 

Let us suppose that lim^ioo j(x) = j±oo ^ 0, then K± 00 = k± OD /j± oa , so if k 2 (x) has 
nice asymptotic behaviour allowing one to define a scattering problem, then so does 
K 2 (x). 

Another possibly more useful substitution (based on what we saw with the Schwarzian 
derivative) is to set J(x)~ 2 = X'(x) with J(x) > 0. We can then write 

K 2 = rlk 2 + ^\. (18) 



k2 = t, \ k2 - i o-r + ? J ir }■ (17) 



1 We should warn the reader that in about 3% of cases "Schwarzian" is mis-spelled as "Schwartzian" . We 
have been guilty of this error ourselves in the past. Worse, this spelling error then sometimes propagates 
and leads to false "folk history". The Schwarz in question is Karl Hermann Amadeus Schwarz (1843 - 1921), 
as in "CauchySchwarz inequality" , with the name of this inequality being incorrectly spelled in some 25% of 
cases. 
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Let us suppose that lim^-too J(x) = J±oo ^ 0, then K± OQ = k± 00 J± 00 , so if k 2 (x) has 
nice asymptotic behaviour allowing one to define a scattering problem, so does K 2 (x). 

These observations about the behaviour at spatial infinity lead immediately and naturally 
to the two key results: 

Theorem 2. Suppose j = 1, (equivalently, J±oc — I)- Then the "potentials" k 2 (x) and 
K 2 (X) have the same reflection and transmission amplitudes, and same reflection and trans- 
mission probabilities. 

Proof: This is automatic since K± OQ = k± OQ , so (jTJ) and the transformed equation f[T3j) both 
have the same asymptotic plane-wave solutions. Furthermore the Miller-Good transforma- 
tion (Q maps any linear combination of solutions of ([1]) into the same linear combination of 
solutions of the transformed equation ({TBI . QED. 

Theorem 3. Suppose j± OQ 7^ \, (equivalently, J± OQ ^ 1). What is the relation between the 
reflection and transmission amplitudes, and reflection and transmission probabilities of the 
two "potentials" k 2 (x) and K 2 (X)? This is also trivial - the "potentials" k 2 (x) and K 2 (X) 
have the same reflection and transmission amplitudes, and same reflection and transmission 
probabilities. 

Proof: The only thing that now changes is that the properly normalized asymptotic states 
are distinct 

explikpox) expliK^x) 

rr — ^ fT7 — ' \ ) 

but map into each other under the Miller-Good transformation. QED. 

3 Improved general bounds 

We already know (from theorem 1) that: 

T > sech 2 \f $dx\ . (20) 



Here T is the transmission probability, and d is the function 

\/(h') 2 + (k 2 -h 2 ) 2 

2ft ' (21) 

with h{x) > 0. But since the scattering problems defined by k(x) and K(X) have the same 
transmission probabilities, we also have 



oc 



00 



T > sech 2 <^ / MX \ , (22) 
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with 
and 



dX = X'dx = jdx, 

^(h x ) 2 + (K*-h 2 y 



2h 



We now unwrap the definitions to yield 

1 \h') 2 + 




2hj 



3 \ 2 j 4 f 



jh 2 



(23) 
(24) 

(25) 
(26) 



This permits us to derive the first form of our improved bound. 
Improved bound 1. V7i(x) > 0,Vj(x) > we have 

"OO 1 



T > sech 



2h 



.3 I 2 j 4 j 2 



dx 



(27) 



Since this new bound contains two freely specifiable functions, it is definitely stronger than 
the result we started from ([2]). The result is perhaps a little more manageable if we work in 
terms of the function J instead of j. We follow the previous logic but now set 



and 



dX = X'dx = J~ 2 dx, 



{h x f + {K 2 - h 2 f 



(28) 



2h 



1 

2h 



+ 



(29) 



This now leads to the second form of the improved bound. 
Improved bound 2. V7i(x) > 0,VJ(x) > we have 

"OO 1 



T > sech 



2h 



(/»') 



J 2 lk 2 



J" 



h 2 
T 2 



dx 



(30) 



A useful further modification is to substitute h = HJ 2 , leading to the third form of the 
improved bound. 
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Improved bound 3. VH(x) > 0,\/J(x) > we have 
T > sech 2 



r°° 1 / 




2 




H' + 2H— 


+ 




J 



J" 

k 2 + — -H 2 

J 



1 2 



dx 



(31) 



Since J{x) > we can without any loss of generality write 



J(x) = exp 



X{x)dx 



(32) 



where x( x ) is rea l but otherwise unconstrained. This permits is to write a fourth form of 
the improved bound. 



(33) 



Improved bound 4. \/H(x) > 0,Vx(x) we have 

T > sech 2 | J ^ sj [H> + 2H X } 2 + [k 2 + X 2 + X' ~ H*] 2 dx} . 

Equations (|27|) . (!30|) . ( 131]) . and (|33|) . are completely equivalent versions of our new bound. 
This improved bound is definitely more general than the original bounds reported in [TTJ 
H~2] . Furthermore, using totally different techniques involving Shabat-Zakharov systems of 
ODEs, we have verified the correctness of these improved bounds by obtaining a completely 
independent proof in references [211 



4 WKB-like bound 

To manipulate this general bound into a form that is "as close as possible" to the standard 
WKB-result, let us first take version 4 of our improved bound, (133]) . and use the triangle 
inequality to write it in a slightly weakened form. 



Improved bound 5. VH(x) > 0, Vx(x) we have 



T > sech 



H' 
2H 



+ X 



+ 



\e + x 2 + x'-H 2 

2H 



dx 



(34) 



We shall choose the functions H(x) and xi x ) to optimize the bound and render it into WKB- 
like form. To do this, introduce a positive parameter A and, similarly to what was done in 
obtaining ((7|), choose H 2 = max{/c 2 , A 2 }. We have in mind here the situation where k 2 (x) 
has a single minimum, and that fc^ in < A 2 < fc±oo, (but still permitting fc 2 iri < 0, so we are 
allowing for the possibility of a classically forbidden region). Then 



T > sech^ 



fc 2 >A 2 



2k 



+ X 



+ 



\x 2 + x'\ 

2H 



+ 



fc 2 <A 2 



1x1 + 



\k 2 + X 2 + x'\ ~ A 2 
2A 



dx 
dx 



(35) 
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Note that if x ~ > 0, then this reduces exactly to ([7]), so the above is a genuine generalization 
of previous results. Now define k 2 = max{0, — k 2 }, so that k — \k\in the classically forbidden 
region k 2 < 0, while k = in the classically allowed region k 2 > 0. 

Furthermore, choose \ = k. Then one has split the integral into three separate regions: 



T > sech 



+ 



fc 2 >A 2 
K + - 



k' 



dx + 



fc 2 <0 . 



2k t 
k'-A 2 
2A 



0<fc 2 <A 2 

dx 



\k 2 -A< 
2A - 



(36) 



The first of these integrals can be explicitly carried out, while the last can again be bounded 
by the triangle inequality This leads to 



T > sech < In 



\k 2 -A 2 



ndx + 

k 2 <0 Jk 2 <0 . 



0<fc 2 <A 2 

W - A 5 



2A 



2A 



dx 



-dx 



(37) 



The last integral appearing above can now be explicitly carried out. Collecting terms and 
rearranging we obtain our WKB-like bound. 



Improved bound 6. Let k^ > AO, then for a single hump potential with the notation 
discussed above 



T > sech 



i i I koo \ ^max AZ/ 

ndx + In — — H — , 

fe2<0 Ui A 2 



0<fc 2 <A 2 



\k 2 - A 2 
2A 



(38) 



Here L denotes the width of the classically forbidden region, and K max denotes the maximum 
height of the barrier. The first term in this bound involves the standard WKB integral J ndx 
over the forbidden region, justifying the appellation "WKB-like". Note that this bound is 
genuinely an improvement over all other published bounds we know of, the closest analogue 
to the current result appearing in our article [1] . That result corresponds to taking A — > koo 
in which case the current bound (|38l) reduces to the less restrictive statement 



T > sech z 



7 . ^max . 
ndx H : h 



kooL 



k 2 <0 



kn 



+ 



\k 2 



k 2 



fc 2 >0 



2k 'r 



-dx 



(39) 



Another key point is that this WKB-like bound ( )38|) contains contributions from the classi- 
cally allowed region. (As in general there must be, potentials with no classically forbidden 
region still generically have nontrivial scattering.) Compare this with the standard WKB 
estimate: 

2 



sech 



ndx + In 2 



fc 2 <0 



(40) 



This or related forms of the WKB approximation for barrier penetration are derived, for 
instance, in many standard textbooks [2j[3llH[5l|6l[7J|8l[9l Under the usual conditions 
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applying to the WKB approximation for barrier penetration we have J ndx 3> 1, in which 
case one obtains the more well-known version 

T WK B ~ exp <^ -2 / Kdx \ . (41) 

I Jk 2 <0 J 

The bound in f l3"51) . and it specialization in f l3"9~j) . are the closest we have so far been able to get 
to obtaining a rigorous bound that somewhat resembles the standard WKB estimate. Note 
the long chain of inequalities leading to these results - this suggests that the final inequality 
( 138|) might not be optimal. We are continuing to search for improvements on this WKB-like 
bound. 



5 Schwarzian bound 



In counterpoint to the WKB-like bound obtained above, there is also a rather elegant bound 
that can in certain circumstances be obtained in terms of the Schwarzian. First, take h = 
(constant) in (I30p . Then 



T > sech z <^ - 
2 



dx 



(42) 



In order for this bound to convey nontrivial information we need lim 2 



>±oo 



JU 2 = h\ oth- 



erwise the integral diverges and the bound trivializes to T > sech (oo) = 0. The further 
specialization of this result reported in [TTJ [T2] and (0J above correspond to taking J = (con- 
stant) = a/^/^oo, which clearly is a weaker bound than that reported here. In the present 
situation we can without loss of generality set h — > k^ in which case 



T > sech 



J2_ 
h 



k 



T 



k 

J 2 



dx 



(43) 



We now need lim^ioo J = 1 in order to make the integral converge. If k 2 > 0, so that there 
is no classically forbidden region, then we can choose J = \fk^fk, in which case 



T > sech z <^ - 
2 



(44) 



This is a particularly elegant bound in terms of the Schwarzian derivative, (|T6|) . which 
however unfortunately fails if there is a classically forbidden region. This bound is also 
computationally awkward to evaluate for specific potentials. Furthermore, in the current 
context there does not seem to be any efficient or especially edifying way of choosing J{x) 
in the forbidden region, so while the bound in fT43|) is both elegant and explicit it is not 
particularly useful. 
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6 Summary and discussion 



The fundamental bounds presented in this note are generally not "WKB-like" - only for one 
particular special case has it proved useful to separate the region of integration into classi- 
cally allowed and classically forbidden regions. In fact, in this WKB-like case it still seems 
that the bound we have derived is sub-optimal, and it is far from clear how closely these 
bounds might ultimately be related to WKB estimates of the transmission probabilities. 
This is an issue to which we hope to return in the future. 

We should mention that if one works with the formally equivalent problem of a parametric 
oscillator in the time domain then the relevant differential equation is 



u(t) + k(t) 2 u(t) = 0, 



(45) 



and instead of asking questions about transmission amplitudes and probabilities one is nat- 
urally driven to ask formally equivalent questions about Bogoliubov coefficients and particle 
production. The key translation step is to realize that there is an equivalence [TJ [TTJ IT2"| 122]: 



1 » T 1-T 
T ^ - -; N <-> 



(46) 



1 + JV' T 

This leads to bounds on the number of particles produced that are of the form 

N < sinh 2 {some appropriate integral} . (47) 

To be more explicit about this, recall that our new improved bound can be written in any 
of four equivalent forms: 



1. For all H(x) > 0, for all x{ x ) we have 



T > sech 



2H 



[H 1 + 2H X ] + [k 2 + X 2 + X'~ H 2 ]dx 



2. For all H(x) > 0, for all J{x) > 0, 

"OO 



T > sech^ 



1 
2H 



H' + 2H 



J>_ 
1 



+ 



T" 

k 2 + — -H 2 

J 



l 2 



dx 



(48) 



(49) 



3. For all h(x) > 0, for all J(x) > 0, 

"OO 1 



T > sech 



2h 



{h'f + 



j 2 U 



T 



j 2 



dx > . 



(50) 



4. For all h(x) > 0, for all j(x) > 

"OO 



T > sech 



2h 



3 \ 2 j 4 f < 



{h>) 



dx 



(51) 
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The equivalent statements about particle production are: 



1. For all H(t) > 0, for all x{t) we have 

N < sinh 2 I.J ^ sj \H' + 2H X f + [k 2 + X 2 + X' - H 2 fdt 

2. For all H(t) > 0, for all J(t) > 0, 



N < sinh 2 



1 

2H 



H' + 2H 



+ 



J" 

k 2 + — -H 2 



dt 



3. For all h(t) > 0, for all J(t) > 0, 



N < sinh 
4. For all h(t) > 0, for all j(t) > 

OO J 

A' < sinlr ( I — 



OO 



J 2 u 2 + 



J" 
T 



/i 2 
J 2 



, 2 _ If 

2 j 4 j 2 



-jh 2 



dt 



(52) 



(53) 



(54) 



(55) 



In closing, we reiterate that these general bounds reported in (12"7}L ( 130|) . ( 13~1T) . and ( 1331 . 
their specializations in (|38|) . ( )39|) . (H3|) . and P4|) . and the equivalent particle production 
bounds in (153|) - (!55|) . are all general purpose tools that are applicable to a wide variety of 
physical situations [131 HH HHJ HS1 El EH]- Furthermore we strongly suspect that further 
generalizations of these bounds are still possible. 
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